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ABSTRACT 

We conduct numerical simulations of axisymmetrical jets expanding into a 
spherical AGB slow wind. The three-dimensional flow is simulated with an axially 
symmetric numerical code. We concentrate on jets that are active for a relatively 
short time. Our results strengthen other studies that show that jets can account 
for many morphological features observed in planetary nebulae (PNs). Our main 
results are as follows. (1) With a single jet's launching episode we can reproduce 
a lobe structure having a 'front-lobe', i.e., a small bulge on the front of the main 
lobe, such as that in the PN Mz 3. (2) In some runs dense clumps are formed 
along the symmetry axis, such as those observed in the pre-PN Ml-92. (3) The 
mass loss history of the slow wind has a profound influence on the PN structure. 
(4) A dense expanding torus (ring; disk) is formed in most of our runs. The 
torus is formed from the inflated lobes, and not from a separate equatorial mass 
loss episode. (5) The torus and lobes are formed at the same time and from 
the same mass loss rate episode. However, when the slow wind density is steep 
enough, the ratio of the distance divided by the radial velocity is larger for regions 
closer to the equatorial plane than for regions closer to the symmetry axis. (6) 
With the short jet-active phase a linear relation between distance and expansion 
velocity is obtained in many cases. (7) Regions at the front of the lobe are moving 
sufficiently fast to excite some visible emission lines. 



1. INTRODUCTION 

The understanding of planetary nebulae (PNs) shaping is the prime focus of the PN 
community for about two decades (Balick & Frank 2002 and references therein). Some 
popular mechanisms for shaping involve jets that are launched by the central star or its 
companion (e.g., Morris 1987, 1990; Soker 1990, 1992, 1996; Soker & Livio 1994; Sahai & 
Trauger 1998; Soker & Rappaport 2000; see Soker & Bisker 2006 for many more references 
and a discussion of the development of the jet shaping models). The shaping jets in these 
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models are blown during the late AGB phase or early post-AGB phase, because mass must 
be supplied to a companion outside the envelope (Morris 1987, 1990; Soker & Rappaport 
2000), or in a post-common envelope evolution (Soker & Livio 1994). Even jets with small 
amount of mass that might form small bullets (termed ansae or FLIERS) are formed in the 
pre-PN phase (Soker 1990). Recent studies strongly support models where the main shaping 
occurs in the pre-PN phase or somewhat earlier (e.g., Sahai et al. 2007). 

Based on many other astrophysical objects, it can be safely deduced that the jets are 
blown by accretion disks (Soker & Livio 1994). A single AGB star does not possess enough 
angular momentum to form an accretion disk, and for that the binary mechanism was in- 
troduced (Morris 1987, 1990; Soker & Livio 1994). In the binary model for the formation 
of bipolar PNs— those with large lobes and with a waist between the lobes— a mass transfer 
process from the AGB or post-AGB star to a compact companion forms an accretion disk 
around the companion (Soker 2002). The companion then blows the two opposite jets. 

Hydrodynamical numerical simulations show that jets launched by the central star (or 
its companion) can indeed account for many features of asymmetrical PNs (Soker 1990; Lee 
& Sahai 2003, 2004; Garcia- Arredondo & Prank 2004; Velazquez 2004, 2007; Riera et al. 
2005; Raga et al. 2007; Stute & Sahai 2007; Akashi 2008; Akashi et al. 2008; Akashi & 

Soker 2008; Dennis et al. 2008; Guerrero et al. 2008). Our main new addition is the usage of 
jets with wide opening angle (Akashi 2008; Akashi & Soker 2008; Akashi et al. 2008), that 
might be more appropriately called CFW, for collimated fast wind. In the present paper we 
explore more properties of lobes inflated by wide jets. We note that rapidly processing jets 
have the same effect as wide jets (Sternberg & Soker 2007), but for numerical reasons we 
simulate wide jets. We simulate cases where the jets is active for a relatively short time, with 
the goal of reproducing the basic type of structures represented by the PN MZ 3 (PN G331.7- 
01.0; 'The Ant') and pre-PN Ml-92 (PN G064.09+04.26; IRAS 19343+2926; 'Minkowski's 
footprint'). The short jet-activity phase leads to a more or less ballistic evolution after the 
jet is shut off, that forms a general linear relation between the distance of gas parcel from the 
center and its outward velocity, as is observed in some PNs (e.g., Mz 3, Santander-Garcia et 
al. 2004; Ml-92, Alcolea et al. 2007; NGC 6302, Meaburn et al. 2008). (This linear relation 
is wrongly termed "Hubble type ffow" . However, in the Hubble flow the space itself expands 
with no center. In a PN there is a well defined center of the expanding gas; see Trimble et 
al. 2007.) 

Because we do not include the ionizing radiation from the central star, we compare our 
results with young PNs, or pre-PNs: The dynamical age of the pre-PN Ml-92 is ~ 1200 yr 
(Alcolea et al. 2008), and the age of the lobes of Mz 3 is ~ 900 yr (Guerrero et al. 2004; 
using a distance of 1.4 kpc from Ghesneau et al. 2007). The young age of the lobes of these 
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objects is one of the reasons we use them for comparison with our results. 



2. THE NUMERICAL SCHEME 

The simulations were performed using Virginia Hydrodynamics-I (VH-1), a high reso- 
lution multidimensional astrophysical hydrodynamics code developed by John Blondin and 
co-workers (Blondin et al. 1990; Stevens et al., 1992; Blondin 1994). We have added ra- 
diative coohng to the code at all temperatures T > 100 K. Radiative coohng is carefully 

treated near contact discontinuities, to prevent large temperature gradients from causing 
unphysical results. The cooling function A(T) for solar abundances that we use was taken 
from Sutherland & Dopita (1993; their table 6). 

We simulate three-dimensional axisymmetrical morphologies. This allows us to use two- 
dimensional axisymmetrical grid, and to simulate one quarter of the meridional plane. There 
are 208 grid points in the azimuthal {9) direction of this one quarter with an equal azimuthal 
spacing, and 208 grid points in the radial direction. The radial size of the grid points increases 
with radius by a factor of 1.015 from one cell to the next. The cells' radial size at the inner 
boundary is Ar(rin) = 2.8 x 10^"^ cm, while it is Ar(rout) = 6.1 x 10^^ cm at the outer 
boundary. In these simulations the grid extends from Tin = 10^^ cm to rout = 4 x 10^^ cm. 

Before the CFW (jet) is launched at t = the grid is filled with a spherically-symmetric 
slow wind having a uniform radial velocity of v{r, 0) = vi = 10 km s"^ The density at t = 
is taken to be 



Ml 



< r < i? 



where i? is a radius that separates two regions: an inner one where the density varies as r~^, 
and an outer one where the density decreases as . The values of (3 for the different runs 
are shown in Table 1. In all new runs presented in this paper we take i? = 3 x 10^^ cm and 
Ml = IO^^Mq yr~^. This density profile represents an AGB mass loss history where the 
mass loss rate increased until a time R/vi before the onset of the CFW (jet), and stayed 
constant at a value of Mi thereafter until t — Q. 

The CFW (jet) is launched from the first 20 radial cell-rows attached to the inner 
boundary of the grid, extending in the range 10^^ cm < r < 6.5 x 10^^ cm, and within an 
angle (half opening angle) a < 9 < a). In the present paper we take a — 30° and CFW 
initial radial velocity of V2 — 600 km s~^ in all runs. For numerical reasons a weak slow wind 
is injected in the sector a <9 < 90°. Here ^ = is along the symmetry axis (vertical in the 
figures) . 
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The initial temperature of the slow wind in this work is Ti = 1000 K. The temperature 
of the gas injected into the CFW is T2 = 1000 K, but it has no influence on the results 
because of the high mach number. The initial temperature of the slow wind has a noticeable 
influences on the final shape of the PN (Akashi 2008). An initial temperature of Ti = 1000 K 
is appropriate for most cases. An initial temperature of Ti = 10^ K might be appropriate for 
some cases that the companion to the AGB (or post-AGB) star is a WD. If the companion 
WD accretes at a high rate it can sustain a constant nuclear burning and becomes super-soft 
X-ray sources (e.g., Starrfield et al. 2004), and ionize the nebula. This is the case for some 
symbiotic nebulae. 

Another important parameter is the duration of the jet active phase. The CFW (jet) 
is blown with unchanged parameters for a time of At2, and then turned off. The main 
motivation for the short duration of the jet active phase is the observed linear relation of 
velocity with distance from the center (e.g., Santander-Garcia et al. 2004; Meaburn et al. 
2008). 

Our calculations do not include the ionizing radiation and fast wind blown by the central 
star during the PN phase. We are studying the shaping that occur in the pre-PN (proto-PN) 
phase. Both the fast wind and the ionizing radiation play a role in the later shaping of the 
nebula (Perinotto et al. 2004; Schonberner et al. 2004), but the basic shaping occurs already 
in the pre-PN or in the late AGB phase (e.g., Sahai et al. 2007). 

3. THE PHYSICAL PARAMETERS 

3.1. The structures to be reproduced 

Our goal is not to explain all features of a particular PN, but rather to show that we 
can reproduce two typical types of bipolar structures as represented by two PNs. The first 
Structure is the one represented by Mz 3 (PN G331. 7-01.0; 'The Ant'), a well studied PN 
(e.g., Guerrero et al. 2004). Although Mz 3 is commonly referred to as a PN, it might 
actually be a symbiotic nebula (Schmeja & Kimeswenger 2001). For our purpose it does not 
matter, as we are looking into the general structure which is shared by many other PNs. 
The image of this PN is given in Fig. [H 

Santander-Garcia et al. (2004) identify four major parts in the nebula. (1) A pair of 
bright bipolar lobes. They extend to a projected distance of ~ 17" (5 x 10^^ cm at a PN 
distance of D = 2 kpc) from the center. (2) Two opposite highly collimated column-shaped 
outflows. These are narrow features along the symmetry axis that are outside the bipolar 
lobes. (3) A conical system of radial structures, also outside the bipolar lobes. (4) A very 



Fig. 1. — HST image of the Mz 3 planetary nebula (Credit: R. Saliai, B. Balick, Hubble 
Heritage Team, ESA, NASA). We mark several morphological features that we will compare 
with numerical results. 



- 6 - 



dim, low-latitude and flattened, ring-like, radial outflow in the equatorial plane outside the 
bipolar lobes. We will try to reproduce only the bipolar lobes, and not the structural features 
outside them, that are the result of previous mass loss episodes. 

There are two other structural components inside the bipolar lobes. A silicate disk of 
radius 0.006" at the center of Mz 3 (Chesneau et al. 2007), and X-ray emitting region near 
the symmetry axis, that might be related to a jet (Kastner et al.2003). The X-ray emission 
and morphology suggest a collimated outflow with a speed of several x 100 km s^^. 

The second kind of bipolar structure we are aiming to produce is the one represented 
by the bipolar pre-PN Ml-92 (PN G064.09+04.26; Minkowski's footprint). An image of 
this PN is given in Fig. [2l Early studies of this pre-PN (proto-PN) were conducted by 
Herbig (1975), Calvet & Cohen (1978), and Trammell & Goodrich (1996), and it was then 
extensively studied by Bujarrabal et al. (1997, 1998a,b) and Alcolea et al. (2007, 2008). 
Because it is a pre-PN, most of the mass is not observed in the visible band. The bipolar 
structure seen in figure Fig. [2] is the illuminated part only; most of the mass resides in an 
expanding molecular disk in the equatorial plane. The visible lobes extend to a distance of 
6" (2.3 X 10^'^ cm at a PN distance of D = 2.5 kpc) from the center. There is a dense central 
disk (not an accretion disk) with a diameter of 10^'' cm, and the total mass of the molecular 
envelope is ~ IM© and the estimated age is ~ 1000 yr (Bujarrabal et al. 1998a,b; Alcolea 
et al. 2008). 



3.2. Previous runs 

Before turning to describe the new runs, we present one of the previously simulated 
cases (Akashi 2008), that also reproduce the basic structure of Mz 3. In that paper different 
mass loss rate histories of the AGB progenitor were studied. In a model presented here in 
Fig. [3l the slow wind mass loss rate was increased by a factor of fci = 6 for a time period of 
Ati = 950 yr, ending at the beginning of the jet-launching phase. We marked on that figure 
morphological features that can be identified in the image of Mz 3 (Fig. [T]). Akashi (2008) 
found that the mass loss history of the AGB wind plays a significant role in the shaping 
process. In the present study the mass loss history is modelled by the density profile of the 
slow wind as given by equation ([1]). 
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Fig. 2.— HST image of the Ml-92 pre-planetary nebula (Credit: ESA, NASA-HST; see 
Bujarrabal et al. 1998b, Trammell & Goodrich 1996, and the Planetary Nebula Image 
Catalogue maintained by B. Balick). We mark several morphological features that we will 
compare with numerical results. 
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3.3. New parameters 

Aiming at the structures represented by the two PNs described in section 13.11 we sim- 
ulated many cases, out of which we present those that are summarized in Table 1. The first 
column lists the run, the second column gives the value of the mass loss rate of one jet M2, 
the third column lists the parameter (3 that appear in equation ([T]), and in the fourth column 
the duration of the jet's launching (active) phase is given. The following physical parameters 
were held the same for all runs: The initial velocity of the slow wind f 1 = 10 km s~^ and 
of the CFW (jet) V2 = 600 km s~^, both of them radial, the initial temperature of the slow 
wind and of the CFW 1000 K, the radius and mass loss rate appearing in equation ([T]), 
R = 3 X 10^^ cm, and Mi = 10~^Mq yr~^, respectively, and the half opening angle of the 
jet a = 30°. 



Table 1: Cases Simulated 



Run 


M2 
{Me yr-i) 




At2 

( yr) 


Figure 


M50 


10-^ 


2 


50 




M3-50 


3 X 10"^ 


2 


50 


m 


M80 


10-^ 


2 


80 




M150 


10-^ 


2 


150 


E 


M600 


10^^ 


2 


600 




M2.2 


10-^ 


2.2 


150 


m 


M3 


10-^ 


3 


150 


[iii[i3i[ii[ia 


M2.2-3 


3 X 10-^ 


2.2 


150 


m 



The first column gives the name of the run, the second column lists the mass loss rate of one side of the 
CFW (one jet), /3 is the density-profile parameter in equation ([1]), and the 4th column gives the duration of 
the CFW injection (jet launching) period. In the last column we list the figures where results of the runs are 
described. In all runs the slow wind velocity is vi — 10 km s^^, the mass loss rate and radius that appear 
in equation ([T]) are Mi = IO^^Mq yr~^ and i? = 3 x 10^^ cm, respectively, the initial temperature of the 
slow wind and the injected CFW is 1000 K, the initial velocity of the CFW is V2 = 600 km s~^, and the half 
opening angle of the CFW is a = 30° . 
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4. RESULTS AND DISCUSSION 

4.1. Dependance on the jet active phase duration 

We run several cases with varying jet-active period At2- One of the motivations for a 
short jet activity period comes from the balhstic motion observe in some PNs, i.e., a hnear 
relation between the distance form the center and the radial velocity (e.g., Santander et al. 
2004; Meaburn et al. 2008). 

We first compare two runs differing only by the mass loss rate of the jet M2. These are 
run M50 presented in Fig. H] and run M3-50 presented in Fig. [5l given each at two times as 
mentioned in the respective figure captions. Run M3-50 evolves faster as more energy was 
injected via the jet. This type of interaction has some instabilities that develop into 'ears' 
protruding from the hot bubble. We can tell that the general structure of the two nebulae is 
the same, but that in run M3-50 the 'ears' are more pronounced. Also, the higher energy in 
run M3-50 manages to prevent the concentration of large amount of gas along the symmetry 
axis. That the two runs have different amount of mass along the symmetry axis suggests that 
this mass concentration is not a numerical effect. It is real, but its exact structure depends 
on the numerics, e.g., the axisymmetric grid forces it to be exactly on the axis. In run M50 
(Fig. H]) the expansion velocity at t = 1140 yr and at a distance of r = 2.7 x 10^'' cm on the 
symmetry axis is Vexp = 58 km s~^, while in run M3-50 (Fig. [5]) the expansion velocity at 
t = 840 yr and at r = 3 x 10^'' cm on the symmetry axis is Vg^p = 85 km s~^. 

These two runs demonstrate how a simple variation in the injected energy can change 
not only the expansion rate, but also some morphological features. This sensitivity is further 
demonstrated by run M80 that is presented at t = 1140 yr in Fig. [61 It is similar to run 
M50, but the jet was active for a time period of At2 = 80 yr instead of At2 = 50 yr in run 
M50. Noticeable is the difference in the hot bubble size obtained in runs M50 and M80. On 
Fig. [H] we try to identify some morphological features appearing in the image of the pre-PN 
Ml-92 (Fig. [2]). There are some quantitative differences in these features, but the general 
structure of the simulations is the same as that of Ml-92, beside the 'ears' that envelope 
the lobes in the simulations but are not observed in Ml-92. Later on we will obtain lobes 
with smooth boundaries, but then other features appearing in Ml-92 are less pronounced, 
or don't exist at all. This discrepancy will have to be resolved when we better understand 
the interaction of the jets with the slow wind. It is possible that a different AGB mass loss 
history will remove the ears if the run is continued to later times. In a future paper we will 
also check whether a more massive jet can form better agreement. 

In any case, we suggest that the basic structure of Ml-92 can be formed by a jet similar 
to that in run M80. The region of the knot along the symmetry axis (or knots if there are 
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several knots) closest to the center of the nebula is within the main shell. This is best seen 
in Fig. [3] (Akashi 2008). This is a result of the usage of a wide jet that form vortices in 
the hot bubble (Fig. [3]). When the jet is narrow, no such knots are formed (Lee & Sahai 
2003). The inner regions of knots cool to low temperatures, but their outskirts can become 
warm, ~ 10^ K (Fig. [7]), and be a source of emission lines, such as found by Bujarrabal et 
al. (1998b) in the knots of Ml-92 (For more on line emission from shocked bullets see Raga 
et al. 2008). 

We cannot predict at this point the X-ray emission of Ml-92. Examining the jet sim- 
ulated here, we find that the shocked jet's material had time to cool (both radiative and 
adiabatic cooling are included in our code) to a temperature of T <C 10^ K. This is seen in 
the temperature map given in Fig. [71 This implies that the shaping jet will not be a source 
of X-ray emission. However, a second jet-launching episode, or a spherical wind from the 
central star, can fill the cavity with a hot enough gas to be a strong source of X-ray emission. 

To close the study of the sensitivity to the duration of the jet active phase, we show 
two more cases in Figs. [SlfTUl In Fig. [H]we present run M150, where beside a vortex around 
the lower density region of the hot bubble, the nebular gas approaches a ballistic motion. In 
Fig. [9] the jet duration is relatively long, and a large bubble with a complicated flow inside 
it is formed. 

The two lobes of Mz 3 might have been formed by two jet-launching episodes (in addition 
to the mass loss episodes that formed the structures outside the lobes), and not by a single 
episode as we examine here. There are some arguments suggesting that two shaping episodes 
of the lobes and front lobes is not unlikely. There have been episodic bipolar ejections before 
the one or two episodes that formed the lobes and front-lobes (Guerrero et al. 2004). The 
X-ray image of Mz 3 suggests the presence of a relatively narrow jet inside the main lobes, 
fitting the size of the front lobe (Kastner et al. 2003). The long-slit spectrum of Mz 3 
reported by Santander-Garcia et al. (2004) shows that the dense partition ring— on the 
boundary of the main lobe and the front lobe— has a broader spectrum. This spectrum 
suggest that this region is moving faster, or that there is a turbulence there. 

In any case, here we point that a single jet-launching episode can form a front-lobe type 
of structure. We marked some morphological features in Fig.[9l that resemble to some degree 
morphological features observed in Mz 3 as marked in Fig. [1] The run presented in Fig. [3] 
(Akashi 2008) better represent the morphology of Mz 3. As stated earlier, our goal is not to 
fit specific PNs, but rather to point to the general morphological features that jets can form. 
In the case of Mz 3 the task is more difficult because the two lobes are not exactly identical, 
in particular in the front lobe region. We see by comparing figures [1] and [9] that the match 
is definitely not perfect, but when the results presented Fig. |3]are considered as well, we can 
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conclude that jets can indeed form the same type of morphological features as observed in 
Mz 3, with different parameters for the two lobes. The lobe appearing in the right side of 
Fig. [T] is better reproduced by the M150 run shown in Fig. [S] (right panel). The 'ears' on 
the front lobe are clearly seen in this numerical run. Also seen is the relatively narrow waist 
in the equatorial plane with a dense ring (torus) at a distance of ~ 6 x 10^^ cm from the 
center. 

As seen in Fig. [HI the flow structure of the gas in the hot low-density bubble, formed 
by the shocked jet material, is quite complicated. For example, we can see a vortex in the 
right side of the bubble, and material flowing inward— a backfiow. The backfiow results from 
the decrease of the pressure in the inner region, which itself result from the expansion of the 
nebula as it follows the dense shell. In Fig. [10] we study the expansion of the dense regions, 
that later form the bright regions of the nebula observed in the visible band. We examine 
only regions having densities of p > 10~^^ g cm~^, and show the magnitude of the velocity 
by a color coding, as indicated in the figure, and the direction by arrows. The front of the 
shell is termed 'cap' (Balick 2004) as marked on Fig. [TOl As evident from the figure, different 
regions in the 'cap' have different velocities that are typically much larger than in the rest 
of the shell. The cap region of Mz 3 also show large velocities (Guerrero et al. 2004). In the 
case of Mz 3 some of the fast moving regions might be outside the main lobe (Guerrero et 
al. 2004), and are not studied here. However, our results show that the 'cap' itself can have 
relatively large velocities. 



4.2. Dependance on the Slow Wind Density Profile 

In Fig. [TT] and Fig. [T^] we present the results for runs where the mass loss rate of the 
slow wind that formed the nebula was not constant. The initial density profile was according 
to equation ([1]) with the value of and other parameters, given in Table 1. 

Two interesting features are seen, (i) Comparing run M3 and M2.2-3 with run M50 
(Fig. H]), we see that the sides of the lobes are smoother, i.e., the 'ears' are less pronounced. 
(ii) A front lobe is seen in both runs M3 and M2.2-3. However, the front lobe in run M2.2-3 
is disrupted, with 'ears' type structure. The type of structure obtained in run M2.2-3 is 
similar in some respects to the disrupted front lobe seen in one of the lobes of Mz 3 (at the 
right side of Fig. [T]). 



4.3. The formation of a torus (ring) 



A dense, slowly expanding ring (torus) is formed in most runs presented in this paper 
and in Akashi (2008). This is an important result proposed by Soker & Rappaport (2000). 
Such rings are observed in many PNs. Huggins (2007) studied several PNs, and found that 
the distance r divided by expansion velocity Vr, t = r/vr, of the jets and the torus in each 
PN give somewhat different values, with lower values for jets. Still, their age is not much 
different and Huggins (2007) concluded that the formation of the jets and the formation 
of the torus are connected. He concluded that the torus (ring) is probably ejected up to 
several hundreds years before the jets. We show that this is not necessarily so. In Fig. [12] 
we show run M3 at two times. We mark several features at t = 540 yr (left panel). We 
then increase the distances of all the marked regions from the center by the same factor, and 
mark the results on the density image at t = 840 yr (right panel), such that the marks on 
the two dense parts along the symmetry axis (the 'jet') coincide. We clearly see that as we 
move from axis toward the equator the dense regions lag behind more and more. Namely, 
it looks as if the jets expands faster than the other regions. The distance from the center 
divided by the difference is smaller for the jet than for other regions, despite that all features 
were formed by the same event. This behavior result from the fast decline in the density, 
/? > 2; performing the same exercise for run M50, where (3 = 2, leads to less conclusive 
results. Namely, the equatorial flow suffers more deceleration than the polar flow, therefore, 
the expansion speed history of different regions is not the same. We did not compare the 
equator itself for a numerical reason. We inject a small amount of material in the equator 
during the jet launching phase for numerical reason. This non-physical material accelerate 
somewhat the ring (torus). 

To further demonstrate this result, we calculate the values of r = r/vr for the regions 
marked in Fig. [13] and few other dense regions, at a later time of t = 960 yr. In Fig. [T4] we 
show the values of r calculated in these regions as function of the angle from the symmetry 
axis. Observational study as that of Huggins (2007) will identify the features at ^ = 0° as 
the jet, and those near 6 = 90° as the torus. It is clear from this figure that the jet appears 
as if it was formed after the torus. However, they were all formed from the same mass loss 
episode. We note that Alcolea et al. (2007, 2008) find the age of Ml-92 to be ~ 1200 yr, 
and claim that it was formed by a single shaping event that last for < 100 yr. 

In a test run (not shown here) we increased the mass loss rate of both the slow wind 
and the jet by a factor of 50. The other parameters were as in run M80. About half of the 
mass of the nebula was concentrated in the torus. We find that we can account for cases 
where the mass in the torus is about equal to the mass in the rest of the nebula. This leaves 
one problem that may still exist with our interpretation of the results of Huggins (2007), 
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and this is, as noted by Huggins (2007), the large mass of the torii in some cases. In cases 
where the torus mass is much more than the mass in the rest of the nebula our model cannot 
account for the whole torus, and a special mass loss episode must be added to the model, 
e.g., by a common envelope interaction. 

4.4. The linear v — r relation of the flow 

As mentioned in section [1] in many PNs there is a linear relation between distance from 
the center and velocity. To check this relation for our simulations, we calculate the radial 
velocity for many small regions in the nebular shell. The radial velocity of each region was 
taken to be the center of mass velocity of the region. In Fig. [15] we draw the radial velocity 
as function of distance for such regions in the M80 and M2.2 runs. The M2.2 run has no 
density map in the paper. It has an initial shallower slow wind density profile than run M3, 
but other than that has all parameters as in run M3. We fit these scattering plots with a 
linear line. A nice linear relation is seen in run M2.2. Although not perfect in run M80, a 
linear relation holds there as well. 

For run M50 and M3 we make the same plots, but mark differently regions according 
to their direction from the symmetry axis, as indicated in the figure caption of Fig. [ini In 
the case of run M50 a general linear relation holds for all regions. In the case of run M3, 
where the slow wind density decline steeply with radius, regions closer to the symmetry axis 
reach higher velocities, as was discussed in section 14. 3[ This is clearly see in Fig. [161 Here 
again the observational results of Huggins (2007) is accounted for by one shaping episode. 
Although there is only one shaping episode, i.e., the jet, this explanation requires that the 
spherical slow wind mass loss rate will increase hundreds of years prior to the jet ejection. 

We conclude that the type of interaction studied here can explain the linear relation 
observed in many PNs, as well as the seemingly result that jets are younger than torii. 

5. SUMMARY 

We conducted numerical simulations of axisymmetrical jets expanding into a spherical 
AGB slow wind. Using the axisymmetrical nature of the flow and its mirror symmetry 
about the equatorial plane, we simulated the 3D flow with a 2D hydrodynamical numerical 
code covering one quarter of the meridional plane. The parameters of the slow wind and 
jets are summarized in Table 1. In the present study we have concentrated on the types of 
structures represented by the bipolar PN Mz 3 and the bipolar pre-PN Ml-92. We used light 
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jets (CFW), i.e., their density is much lower than the slow wind density at the same radius. 
Our goal was not to play with the different parameters in order to fit as best as possible each 
structure, but rather to show that the basic types of structures can be reproduced, and to 
study the role of some parameters. As indicated in the caption of Table 1, many parameters 
were held fixed in the present study. 

We did not vary the slow wind mass loss rate at the end of the AGB, Mi = 10~^Mq yr~^, 
that appears in the density profile (eq. [T]). By scaling the mass loss rate in the slow wind and 
the jet one can produce more massive nebulae. In a test run (not shown here) we increased 
the mass loss rate of of both the slow wind and the jet by a factor of 50, keeping the other 
parameters as in run M80. We obtained a nebula similar in structure to that in run M80, 
but typical densities larger by a factor of ~ 50 in the different regions, and much lower 
temperatures. The main differences in the structure were the structure of the ' they 
are the results of instabilities that are sensitive to other small vatiations. These differences 
are the result of the shorter radiative cooling time of the denser gas. 

Our main results are summarized below. One should remember that the axisymmetrical 
numerical grid forces some features, in particular dense clumps, to be formed exactly along 
the symmetry axis. These are real features, but in reality such features can be off-axis and 
more dispersed. 

1. With a single episode of jet's (CFW) launching we could reproduce a lobe structure 
having a 'front-lobe', as indicated in Figs. [9] and [3l the run represented in Fig. [3] is 
from Akashi (2008), and does not appear in Table 1. the 'front-lobe' is similar to the 
one seen on one of the lobes of Mz 3, as indicated in Fig. [T] (left lobe). The other 
lobe of Mz 3, on right side of Fig. [T], has a structure that was reproduced by the run 
shown in Fig. [8], and to some degree by that shown in Fig. [121 We cannot rule out the 
possibility that the front lobe was formed by a separate jet; we only showed that it 
can be reproduced by a single jet-launching episode. 

2. In some runs (e.g.. Fig. [6]) dense clumps are formed along the symmetry axis, such as 
those observed in the pre-PN Ml-92 (Fig. [2]). 

3. Instabilities lead to the formation of 'ears' on the lobe (bubble) boundary, e.g.. Fig. [61 
Such 'ears' are observed in the front-lobe of one of the two lobes of Mz 3 (on the right 
side of Fig. [T|), but in many other PNs such ears are not seen. 

4. We run several cases where the density of the slow wind declines faster than r~^, i.e., 
/3 > 2 in equation ([T|), as is expected when the mass loss rate increases at the end 
of the AGB. In those cases (Figs. [TH [T2D the 'ears' disappear from most of the lobe 



- 15 - 



boundaries, and might appear only in the front. In particular, in Fig. [12] we can see a 
structure where only the front lobe is disrupted by 'ears'. 

5. The different runs performed here and in Akashi (2008) show that the mass loss history 
of the slow wind has a profound influence on the PN structure. Relatively small 
variations in the slow wind mass loss history can lead to new morphological features. 

6. A dense expanding torus (ring; disk) is formed in most of our runs, as proposed by 
Soker & Rappaport (2000). It is likely that in some types of binary interaction an 
expanding torus is formed as well. Our results show that a ring (torus) is easily formed 
by jet shaping, with no need for a specific equatorial mass loss process. However, if 
the torus is very massive, much more than half the total mass in the nebula, then a 
special equatorial mass ejection process should be added to the model (Huggins 2007). 

7. The torus and lobes are formed at the same time and from the same mass loss rate 
episode. However, in runs with (3 > 2 the ratio of the distance divided by the radial 
velocity, r = r/w^, is larger for regions closer to the equatorial plane (fig. [T3|, [T ^ [T6l) . If 
the ratio r = r/f^ is taken to be the age of the specific regions, then one might deduce 
that the torus was formed before the jets, although this is not the case: they were all 
formed from the same shaping mass loss episode. 

8. Our results, as depicted in Fig. [T^ and [TBI show that impulsive jets, i.e., a short active 
phase, lead to a linear relation (wrongly termed by some authors as a Hubble flow) 
between distance and expansion velocity. 

9. Dense regions in the shell reach sufficiently high velocities to shock-excite visible emis- 
sion lines (Fig. [To]). Indeed, some dense regions reach temperatures of > lO'' K (Fig. [7]). 

We thank Bruce Balick, Patrick Huggins, and Raghvendra Sahai for many clarifying 
and useful comments. This research was supported in part by the Asher Fund for Space 
Research at the Technion. 
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Density map and velocity field 




Fig. 3. — Results from a run conducted by Akashi (2008; it does not appear in Table 1 
here). Shown are the density and velocity map, and the temperature map, both maps at 
t = 1560 yr. In all runs the jet (CFW) starts at t = 0. The initial (at t = 0) slow wind velocity 
is f 1 = 10 km s~^, and the density profile is built for a mass loss rate of Mi = 3 x 1O~^M0 yr~^ 
for r < 3 X 10^^ cm, and Mi = 5 x 1O"''M0 yr"^ for r > 3 x 10^^ cm. The jet has a half 
opening angle a = 30°, it is injected with a radial velocity of V2 = 600 km s~^, and a mass 
loss rate into one jet of M2 = 1O~^M0 yr~^. The initial temperature of the slow wind and 
the injected jet material is 10000 K, and the temperature is limited from below at 1000 K. 
Marked are features that can be identified in the image of Mz 3 (Fig. [1]). Arrows indicate 
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Fig. 4. — Density map for the M50 run at t = 540 yr (left) and t = 1140 yr (right). The time 
is set to if: = when the jet starts. The density scale here and in all figures is in logarithmic 
units of g cm~^. 
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Fig. 5. — Density map for the M3-50 run at t = 540 yr (left) and t = 840 yr (right). 
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Fig. 6. — Density map for the M80 run at t = 1140 yr. Marked are some morphological 
features, some of which can be identified in the image of Ml-92 (Fig. [2]). 
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Fig. 7. — Temperature map of the M80 run at t = 1140 yr. The temperature scale on the 
right is in logarithmic units of degrees in Kelvin. 




Fig. 8. — Density map for the M150 run, at t = 540 yr (left) and 1140 yr (right). 
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Fig. 9. — Density map for the M600 run aX t = 840 yr. Marked are features that can be 
identified in the image of Mz 3 (Fig. [1]). Arrows indicate flow direction, with a scale of: 
V > 200 km s~^ (long arrow), 20 < f < 200 km s~^ (medium arrows), and f < 20 km s~^ 
(short arrows). 
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Fig. 10. — Velocity map for regions having a density of p > 10~^^ g cm~^ for the M600 run 
at t = 840 yr. The magnitude of the velocity is represented by the color, with the scale as 
given by the bar in units of km s~^. The arrows give only the direction. The rectangle show 
the region termed 'cap'. 
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Fig. 11. — Density map for the M3 run at t = 840 yr. 
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Fig. 12. — Density map for the M2.2-3 run at t = 840 yr. Arrows indicate flow direction, 
with a scale of: v > 200 km (long arrow), 20 < f < 200 km s^^ (medium arrows), and 
w < 20 km s~^ (short arrows). 



- 25 - 




Fig. 13. — Density map for the M3 run at t = 540 yr (left) and t = 840 yr (right). Several 
dense regions are marked on the left panel. The distances of the marked regions are increased 
by the same factor, and copied to the right panel to match the location of the jet. The other 
regions 'lag' behind, with differences increasing toward the equator. Namely, the jet is 
moving relatively faster. Dividing the distance from the center by velocity for the different 
regions might give the wrong impression that the jet was launched after regions closer to the 
equatorial plane. 
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Fig. 14. — The values of r = r /v^ for the regions marked in Fig. [T3] and few other dense 
regions at a later time of t = 960 yr. Here Vr is the radial velocity of the region, and r its 
distance from the center. 6 is the direction of the dense region measured from the symmetry 
axis; 6' = is the jet direction. Although all features where formed by the same mass loss 
episode, it appears as if the jet is younger than features formed closer to the equator (at 
e = 90°). 
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Fig. 15. — Radial velocity (circles) and total velocity (dots) versus distance from the center 
of different small regions in the dense nebula for the M80 (left) and M2.2 (right) runs. 
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Fig. 16. — Radial velocity versus distance from the center of different small regions in the 
dense nebula for the M50 (left) and M3 (right) runs. Dots marked region that are within 
25° from the symmetry axis (vertical axis), while asterisks and squares mark regions within 
a direction 25° — 50° and 50° — 75° from the symmetry axis, respectively. 



